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cn : 

■ The Cornwall- Jackiw-Tomboulis (CJT) effective action for composite opera- 

tors at finite temperature is used to investigate the chiral phase transition within 
the framework of the linear sigma model as the low-energy effective model of 
O |' quantum chromodynamics (QCD). A new renormalization prescription for the 

D \ CJT effective action in the Hartree-Fock (HF) approximation is proposed. A nu- 

00 ■ merical study, which incorporates both thermal and quantum effect, shows that 

1^ | in this approximation the phase transition is of first order. However, taking into 

account the higher-loop diagrams contribution the order of phase transition is 
t-H ' unchanged. 

> : 

\jD , PACS numbers: 11.10 Wx, 11.10 Gh, 11.30 Rd, 05.70 Fh. 

m ■ 
o ■ 
o\ : 
o . 

m ■ I. INTRODUCTION 

O ■ 

In recent years the study of the phase transitions has become the subject of intense 
investigation since it is very important from various aspects, ranging from cosmology to 
particle physics. According to the big bang model there is a series of phase transitions, 
^ ; including, of course, the chiral phase one, at the early stage of the universe evolution. 

The nature of the electroweak phase transition has called much attention due to the 
^ . suggestion that the baryon asymmetry may be generated at the electroweak scale if the 

transition is of first order [1]. The restoration of chiral symmetry may lead to several 
interesting phenomena observed, for instance, in the dilepton mass spectrum [2] or 
in the formation of disoriented chiral condensates [3]. The quantum chromodynamics 
(QCD) is chirally symmetric for light quark sector and must be spontaneously broken 
due to the existence of pions, the Nambu-Goldstone (NG) bosons. The lattice QCD 
calculations indicate that the restoration of chiral symmetry occurs at temperatures of 
order 150 MeV [4], which is expected to be probed in ultrarelativistic heavy-ion collision 
experiments planned at RHIC and LHC [5]. The field theoretical investigation of the 
symmetry restoration at finite temperatures was first carried out by Kirzhnits and 
Linde [6] and then was systematically developed by Weinberg [7] , Dolan and Jackiw [8] 
as well as many others adopting the effective action as the most appropriate formalism 
for such studies. It was proved that the leading contributions at very high temperatures 
come from infinite series of certain class of multiloop diagrams in perturbation theory, 
which, in the A0 4 theory, are the so-called daisy and superdaisy diagrams. Recently 
the effective action for composite operators, originally introduced by Cornwall, Jackiw 
and Tomboulis [9] at zero temperature, has been extended for finite temperatures in 
A0 4 theory by Amelino-Camelia and Pi [10], who showed that one needs to resume 
only "double bubble" graphs instead of summing infinite set of daisy and superdaisy 
graphs using the tree level propagators. This is known as Hartree-Fock approximation 
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[11]. Moreover, the CJT formalism at finite temperatures leads to reliable results [12, 
13], among other approximate approachs. 

Many works on chiral symmetry restoration at high temperatures have been accom- 
plished [14-20] within the framework of the linear sigma model, which is considered to 
be the relevant model [21] for effective theory in low-energy phenomenology of QCD 
modelling the hadron dynamics associated to the chiral phase transition. However, 
there exists serious difficulty concerning the renormalization of the CJT effective ac- 
tion in the Hartree-Fock approximation. The point is that a consistent renormalization 
cannot be performed in the broken phase because all terms in the effective action are 
expressed not only by renormalized quantities but bare quantities remain in several 
terms [15]. In [16,19] one tried to regularize divergent integrals by means of several 
schemes instead of renormalizing the CJT effective action in Hartree-Fock approxima- 
tion at finite temperatures. 

In this work, dealing with the chiral phase transition at finite temperatures, we re- 
examine the linear sigma model since it is most suited for our purpose of describing 
the phase transition in a self-consistent approximation. We do not consider the general 
O(N) linear sigma model at large-N limit because it significantly simplifies the counting 
of multiloop bubble diagrams, which is the most common source of inaccuracy. 

We use the imaginary time formalism of Matsubara [22] and work in Euclidean space- 
time. The Feynman rules are the same as those at zero temperature, except that 

where u) n = 2nn/j3. 

This paper is organized as follows. Section II is devoted to calculating the CJT 
effective potential at finite temperatures in the Hartree-Fock approximation. Here a 
new renormalization prescription is proposed. A numerical computation is presented in 
Section III. The higher-loop diagram contribution in the vicinity of critical temperature 
is considered in Section IV. Finally conclusion and discussion are given in Section V. 



II. THE CJT EFFECTIVE POTENTIAL 



The lagrangian density of the 0(4) linear sigma model reads 

C=\(d^f + ^ + ^)\ (11.1) 

where 2 = (f) a (f) a ; ft = n\ % = 1, 2, 3; 4 = a. 
The counterterms, which must be added to (II. 1), are chosen as 

A£ = ^0 2 + ^(0 2 ) 2 + ^, (IL2) 

in which the counterterm 5fl is introduced for vacuum energy. As a rule, all divergencies 
of the theory have to be absorbed by the counterterms (II. 2). 

By shifting the field as 4> a (x) — > 4> a (x) + 0q > the tree-level propagators D a p have the 
form 
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8(j) a {x)5(j)P{y) 



4>=<j>o 



d 2 + m 2 + 5m 2 + ^ ' U 2 
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x 



Here J(0) is the classical action 



= J d 4 x(£ + A£). 
In momentum space, is written as 



k 2 + m 2 + bra 2 4 A - 2 
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5^ + 



(II.3) 



(II.4) 



The interaction action, describing the vertices of the shifted theory, is given by 

A + SX , 9 ,„ ,„ A + tt 2 2 



d 4 x 



i 2 « + 



Then, for constant <po, we arrive at the CJT effective potential at finite temperatures 
V^ JT (cPo,G) = J(0o) + l -Jti\nG-\k) + \Jti [D-'ik^oWk^+Vf^i^G). 

Here G a p(k) is the full propagator of the theory and V^ JT is the sum of all two and 
higher order loop two-particle irreducible vacuum graphs of the theory with vertices 
given by I int ((f), (f> ) and propagators set equal to G a p(k). The graphs shown in Fig. 1 
are under discussion. It is clear that, among them, only the two-loop graph, the "double 
bubble" one, of order 0(A) includes contributions from daisy and superdaisy graphs 
of ordinary perturbation theory. It is known that truncating the series at 0(A) is the 
Hartree-Fock approximation. 








FIG. 1. Two-particle irreducible graphs contributing to Yq.{<P, G) up to the three- loop level in a Xtf) 4 theory. 
The solid line represents the propagator G(x,y). There are two kinds of vertices: a three-particle vertex 
proportional to \<j> and a four-particle vertex. 



Therefore, the CJT effective potential at finite temperatures in the Hartree-Fock 
approximation is derived 

v--(0 o , G) = sn + + X + 6X '^ 2 
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+i J tr In G~\k) + -J tr [D~\k- O )G] 



+ 



'P 
X + 5X 
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G aa {k>) 



+ 2 / G al3 (k) / Gp a (k) }. 
J p J p 



(II.5) 



Minimizing Vf JT (<f>o,G) with respect to G a p(k) we obtain the Schwinger-Dyson (SD) 
equations for propagators 



X + 5X 
~~ 6 



S a p [ G ss (k) + 2 ! G a p(k) 

J /3 J /3 



(II.6) 



Basing on the structure of (II. 6) and the 0(4) symmetry we adopt the following ansatz 
for propagators G a p in Hartree-Fock approximation 



G- 1 p(p) = (p 2 + X")5 a , + Y' 



(II.7) 



Inserting (II. 7) into (II. 5) provides 



Vf JT (<f> , M„ M w ) = 5n+ m2 \ 8m \ l + ^(^) 2 + Q(AV) WW* ) 



where 



24 

+ ; [ m 2 + 5m 2 + ^^<Pl - Ml ) P(M a ) 
+ nr ■ 5m 2 ■ ^-A>f,- M 2 )P 
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X + 5X 
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{3 [P(M CT )f + 15 [P(M u )f + 6P(M ff )P(M w )} , (II.8) 



M a = X 2 + Y 2 , 
M w = X 2 , 



and 



Q{M)= l -J\n{k 2 + M 2 ) ) 
P(M) = 



1 



k 2 + M 2 ' 



In order to regularize the divergent integrals P(M) and Q(M), contained in the 
expression (II. 8) of the effective potential, let us use the three-dimensional momentum 
cutoff scheme, in which every divergent integral is written as the sum of a divergent 
part and a finite part, namely 

Q(M) = DivQ(M) + Qf(M), 



M 4 M 2 
DivQ(M) = - — h + ^h, 



d 3 k , / g(fc)\ 
- — — In 1 - e t 
(2vr) 3 V J ' 



P(M) = DivP(M) + P/(M), 

DivP(M) = h - M 2 I 2 , 

M 2 M 2 f d 3 k 
P/(M) = — In- 



16tt 2 



/// 



P(A;)(l-e^) 
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Now the renormalization is carried out so that all infinite terms appearing in (II. 8) 
must be eliminated. For this end, the counterterms SQ, 5m 2 and 5\ are chosen to obey 
the following constraints: 



Sn + DivQ(M a ) + 3DivQ(M 7r ) = 0, (II.9) 
2 + P(M a ) + 3P(M 7r )] 



<5m 2 



2 



5A / (^0) 2 , j2 r>( Tl/r \ , ^nM, \ , 1 |D ,» N1 2 , 5 rn/ _ Nl2 



+ T \ 6 + ^ P < M ") + ^P(M.) + - [P{M a )Y + - [P{M^)Y + P{M a )P{M 7T ) 
+i (m 2 + ^0 2 - M 2 ) DivP(M CT ) + ^ (m 2 + ^0 2 - M 2 ) DivP(M 7r ) 

{ 3 Div [P(M CT )] 2 + 15Div [P(M n )] 2 + 6Div [P(M ff )P(M w )]} = 0, (11.10) 

(11.10) is a linear equation of two unknowns Sm 2 and 5 A, it has an infinite number of 
roots. In the following are two simplest cases: 

a) 5m 2 = 0, 

b) 5\ = 0. 

In the first case 5X tends to a finite limit as A — > +oo and in the second case 5m 2 
tends to an infinity as A — > +oo. As a rule, both cases are equally accepted for 
renormalizing the CJT effective potential in Hartree-Fock approximation. In general, 
the counterterms are temperature dependent as was mentioned in [14]. 

We obtain ultimately the renormalized effective potential at finite temperatures 



m 2 j2 A 

1 / 2 , A 

m 



2 V 2 



if JT (0 O , M a , M w ) = + _(0 2 ) 2 + Q f (M a ) + SQ^AQ 

$ " M 2 ) P,(M CT ) + | (m 2 + j^0 2 - M 2 ) P,(M.) 
+^ {[P^M.)] 2 + 5 [P/tikQ] 2 + 2P / (M (J )P / (M 7r )} . (11.11) 

The SD equations for M 2 and M 2 and the gap equation for the sigma condensate are, 
respectively, derived from (11.11) for 0q = a, O = n 1 = 0, 

Ml = m 2 + + ±P f (M a ) + ^P f (M n ), (11.12) 

A A ^ A 

M 2 = m 2 + -a 2 + -P f (M a ) + -P f (M„), (11.13) 



and 



V rf 0o / 0q = ff o 

00 = 



or 
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m 2 + + ±P f (M a ) + \P{M„) = 0. (11.14) 
Taking into account (11.12) the equation (11.14) reduces to 

^ 2 -M 2 = 0. (11.15) 

Adopting the argument of [19] we distinguish the effective masses of mesons, M a 
and M n , as those appearing in the approximate propagators, from the physical masses 
which are defined as the curvatures of the effective potential around its minimum, 



^ 2 def. 

da 2 



(11.16) 

a = <7 



7V = 

2 def. fdPV^ 



dn 2 



= 0, (11.17) 
a = <7 

7T = 



in the broken symmetry phase. 

Equations (II. 12)-(II.15) enable us to determine the "transition temperature" T cl , 
at which the broken chiral symmetry is restored: <7o, M CT and M n tend to zero as 
T -> T cl - 0. Indeed, Eq. (11.15) tells that M a -> as tr -»■ and Eq. (11.13) simplifies 
to 

M 2 = m 2 + \Pf(M n ). (11.18) 
When M„. vanishes as T — > T cl , ie. M^/T c i <C 1, we have the expression 



P/(M W )«7^ 



cl 



1 M, 1 M 2 M 2 
In 



12 4tt T c i 16tt 2 T 2 



(11.19) 



Combining (11.18) and (11.19) leads to 

,2\ V2 



T cl « v 7 ^ ( — j = v^/tt « 131.5 MeV. 



Here f w = 93 MeV is the pion decay constant. 



ill. NUMERICAL STUDY IN HARTREE-FOCK APPROXIMATION 

Let us first introduce a term explicitly breaking chiral symmetry into the lagrangian 
density 

C = \ (d^f + + + «J - ~. (III.1) 

The renomalized effective potential corresponding to (III.l) reads 



G 



771^ A 

V p c JT (a , M(j, M w ) = — al + -a 4 - ar + Q/(M CT ) + 3Q f (M n ) 



+\ (- 2 + \°l - P f {M a ) + ? (m 2 + ^ - M 2 ) Pf(M n ) 
+ ^ {[P/(M CT )] 2 + 5 [P/^)] 2 + 2P / (M (7 )P / (M 7r )} , (III.2) 
which leads to the gap equation for the sigma condensate o"o 

a (V 2 - ±o£) = c. (III.3) 

The SD equations for M a and M n , derived from (III. 2), have the same form (11.12) 
and (11.13). Inserting these equations into (III. 2) we arrive at 

77l 2 A 

V(a , T) = —a 2 + -4 - ca + Q f (M a ) + 3Q f (M w ) 

" g [^/(M CT )] 2 - - [P/(M.)] 2 - I P / (M.)P / (M 7r ). (III.4) 

Now the physical mass of pion no longer vanishes in broken symmetry phase and 
equals to 

ml = — . (III.5) 

For numerical computation we use the model parameters at zero temperature as 
initial condition, namely 

c = Uml at T = 



A = ^(m 2 -m 2 ) atT = 
1 o 3 



m 2 = — m 2 + -ml at T = 0. 

We take m w = 138 MeV, m a = 500 MeV and f n = 93 MeV at T = 0. 
For convenience, let us consider two separate cases. 



A. Hartree-Fock approximation in chiral limit c = 

In addition to the preceding parameters, the regularization introduced another pa- 
rameter /i 2 , which corresponds to renormalization scale. Therefore, we must first de- 
termine a suitable value for /i 2 , then all other quantities are determined at jj? = /i 2 ,. 
The value /i 2 , is determined as the real root of the following equation 

a (/x 2 ,0)| M2=M 2 = / 7r = 93 MeV, 

where a (/i 2 ,0) is a solution of the system of equations (11.12), (11.13) and (11.14) at 
T = 0. 
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The numerical computation gives /io = 539.6 MeV 1 . Eliminating o"o from (11.12)- 
(11.14) we get 

Ml = -2m 2 - XP f (M a ) - XP f (M n ), 

Ml = ~P f (M a ) + ^P f (M n ). (III.6) 
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FIG. 2. The temperature dependent solution of 
the SD equations in the chiral limit. 



FIG. 3. Evolution of the order parameter as a 
function of temperature. 



Inserting /x = /x into the system (III. 6) and then solving numerically the preceding 
system of two equations we obtain the solution presented in Fig. 2. As is shown 
in Fig. 3 the temperature T c i, calculated above by means of the high temperature 
approximation, is very close to the one, obtained by numerically solving the system 
(II.12)-(II.14). As the temperature increases from zero, the order parameter decreases 
from f w , jumps to zero at T cl and remains zero above T cl . In the meanwhile, the 
effective masses of mesons, M a and M n , change along the upper lines, jump to lower 
coinciding line at T& and increase along this line. It is clear that the phase transition 
is of first order. The indication of a first order phase transition has been reported in 
many publications [14-20]. 




20 40 60 80 100 

So 

FIG. 4. Evolution of the effective potential V(ao,T) as a function of the order parameter ao for several 
temperature steps: T = 197, 193, 189, 185, 181, 177, 173 MeV from top to bottom. Two minima of V(a ,T 
appear as degenerate at T c ~ 185 MeV. 



x which corresponds to M CT = 481.58 MeV, M„ = 73.53 MeV at T = 0. 
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Next let us compute numerically the effective potential V(ao,T), given by (III. 4), as 
a function of the temperature and the order parameter. This will give more insight into 
the nature of phase transition. In Fig. 4 is depicted the evolution of V(a ,T) against 
o"o for several temperature steps. It is found that the two minima of V(ao, T) appear as 
degenerate at T c2 pa 185 MeV. For T > 193 MeV, V((Tq,T) has only one minimum at 
o"o = 0. The shape of the potential confirms that a first order phase transition occurs. 



B. Hartree-Fock approximation in the broken symmetry case c ^ 



For completeness let us now consider the case c ^ 0, hence the gap equation for the 
order parameter (11.14) turns out to be 



cr - C = 0, 



(III.7) 



and the SD equations (11.12), (11.13) remain unchanged. 

After solving numerically the system (11.12), (11.13) and (III. 7) we get the solution 
presented in Fig. 5. 
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FIG. 5. (a) Solution of the system of gap equations for c / 0. At low temperature M n coincides with the 
observed m n . (b) Evolution of the order parameter as a function of temperature. 



Comparing the above presented result with those obtained in [16-19] we have an 
excellent agreement which shows that the quantum effect is an negligible quantity. 



iv. THE HIGHER-LOOP DIAGRAMS CONTRIBUTION 



In the preceding section it was shown that the phase transition is of first order. The 
Gaussian approximation [24] in 3+1 dimensions gives the same result. However, the 
renormalization group approach applied to the linear sigma model [25, 26] indicates 
that the phase transition is of second order. Arnold and Espinosa [27] pointed out that 
loop diagrams other than superdaisy ones are important near critical temperature. 
The nature of the phase transition in the 4 field theory remains a basic question to 
be settled [28]. Therefore, this requires a further investigation on higher-loop effect. In 
this respect, it is necessary to incorporate the higher- loop diagrams into consideration. 
The higher-loop graph next to the double bubble one is the sunset graph given in Fig 6 
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FIG. 6. The sunset graph. 

In a recent publication [29] Bordag and Skalozub studied, in addition to the double 
bubble, an infinite series of graphs shown in Fig 7 

^fT>rnwn+^frf) 

2 12 V / 48V\ / / ^2" \j ^2 n+l n \ 



FIG. 7. The big dots symbolize further insertion of the subgraph — 0~ so that the number of vertices is n. 
The 3-particle vertex ~ — \4>o, the 4-particle vertex ~ —A. 

which is motivated by the 1/N expansion. 
Making use of the ansatz 

G-\p) = p 2 + M 2 

at high temperature they concluded that the contribution of this set of graphs does 
not influence significantly on the physical process, making the transition a bit stronger 
first order. 

Next let us generalize the consideration including an infinite series of diagrams de- 
picted in Fig 8 




FIG. 8. The big dots symbolize further insertion of the subgraph — 0~ so that the number of vertices is n. 

Then the effective potential, corresponding to Fig 8, reads 
where 

H = Tr 9 G(g)S(g) {h[XE(q)\ + / 2 [AS(g)]} , 
K = Tr g S 2 (g) {#i[A£(g)] + £ 2 [A£(g)]} , 
Z(p) = Tr q G(q)G(p + q) = -0- 

and fi, f2,gi,g2 are respectively the functions representing the series expansions given 
in Fig 8 with the coefficients a n ,b n ,c n ,d n . 
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It is evident that 

f i [0]=g i [0] = l, 1 = 1,2. 

Now it is easily seen that we arrive at the same conclusion by applying the discussion 
of [29] to every fixed expression for every function among f\, f2,gi and g 2 - 

v. CONCLUSION AND DISCUSSION 

We have considered the 0(4) linear sigma model at finite temperature within the 
framework of the CJT formalism restricted to the double bubble diagram approxi- 
mation. The renormalization of the CJT effective potential for this model has been 
discussed in many papers [14-19] and one succeeds only in the large- N limit case [23], 
where all divergent terms are absorbed into the bare quantities. For other case, how- 
ever, not all divergences can be removed away [15]. In the present paper, following 
exactly the spirit of renormalization that requires all divergent terms must be ab- 
sorbed into counterterms, corresponding to renormalizing mass and coupling constant, 
we impose two constraints (II. 9) and (11.10), to ensure that only finite terms would be 
still present in the effective potential. It is possible that this renormalization method 
is most suited for renormalizing the CJT effective potential in loop approximations. 
In order for the Goldstone theorem not to be violated the physical masses of mesons 
are defined by (11.16), (11.17) and propagators are chosen in the symmetric form (II. 7). 
The solution of the gap equations, derived from the renormalized effective potential, 
and the shape of the effective potential, as a function of the order parameter exhibits 
a first order phase transition with the critical temperature T c ~ 185 MeV. 

To go further to higher-loop graphs we have shown that the order of the phase 
transition is unchanged by incorporating an infinite series of diagrams, in which the 
sunset one is a special case. 

In summary, it is possible to confirm that the phase transition in the 4 field theory 
is of first order. 
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